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Abstract. We investigate a driven diffusive lattice gas model with two oppositely 
moving species of particles. The model is motivated by bi-directional traffic of ants 
on a pre-existing trail. A third species, corresponding to pheromones used by the ants 
for communication, is not conserved and mediates interactions between the particles. 
Here we study the spatio-temporal organization of the particles. In the uni-directional 
variant of this model it is known to be determined by the formation and coarsening 
of "loose clusters". For our bi-directional model, we show that the interaction of 
oppositely moving clusters is essential. In the late stages of evolution the cluster 
size oscillates because of a competition between their 'shredding' during encounters 
with oppositely moving counterparts and subsequent "coarsening" during collision- 
free evolution. We also establish a nontrivial dependence of the spatio-temporal 
organization on the system size. 
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1. Introduction 

Systems of interacting driven particles are of current interest in statistical physics 
for understanding the fundamental principles of spatio-temporal organization far from 
equilibrium as well as for possible wide ranging practical applications PUI2]- The simplest 
model of this type with only one species of particles is the so-called totally asymmetric 
simple exclusion process (TASEP). In the TASEP, a particle is picked up randomly and 
moved forward by one lattice spacing, with the hopping probability q + , provided the 
target site on the one-dimensional lattice is empty 

Most of the earlier models with more than one species of particles [3J H EJ El HI El El ED] 
addressed mainly the questions on the (im)possibility of spontaneous symmetry 
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Figure 1. A snapshot of an ant-trail on a hanging cable. It can be regarded as strictly 
one-dimensional. But, nevertheless, traffic flow in opposite directions is possible as two 
ants, which face each other on the upper side of the cable, can exchange their positions 
if one of them, at least temporarily, switches over to the lower side of the cable. 

breaking, phase separation, etc. in one- dimensional stochastic driven systems. Here 
we extend the two-species models like, for example, ABC model jH |3] and AHR model 
jH] to develop a three-species model [TTj. The density of the new species introduced in our 
model, which corresponds to the pheromones secreted by the ants for communication, 
is not conserved by the dynamics of the system and its coupling to the two other 
conserved variables leads to uncommon features of the spatio-temporal organisations of 
the conserved particles. For a special choice of the model parameters our model becomes 
identical to the AHR model classes (see Sec. 3). 

Our model is motivated by a real-life situation where bi-directional ant-traffic in an 
effectively one-dimensional system is observed. Fig. ^ shows a hanging cable which is 
part of an ant-trail. Ants prefer moving on the upper side of the cable, which would 
make the motion strictly one-dimensional. If ants moving in opposite directions meet 
head-on, after short "negotiations" one of them would switch to the lower part of the 
cable temporarily in order to make passing possible. Due to the "negotiations" this 
process leads to slowing down of both the ants. We capture this entire process by 
an exchange of the positions of the two ants with a sufficiently low probability which 
corresponds to a slower movement of each of them as compared to a freely moving ant. 
The mathematical formulation of the model in terms of such exchange probabilities is 
very similar to the formulations of the class of one-dimensional stochastic models of 
driven-diffusive lattice gases to which the AHR and ABC models belong. 

The number of particles leaving a site per unit time is called the flux or current; the 
flux-versus-density relation is often referred to as the fundamental diagram. We study 
interplay of spatio-temporal organization and flow (fundamental diagram) in the system 
by a combination of analytical arguments and extensive computer simulations. 
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Figure 2. Nontrivial transitions and their transition rates. Transitions from initial 
states PL, 0L and OP are not listed. They can be obtained from those for LP, L0 and 
PO, respectively, by replacing R <-> L and, then, taking the mirror image. 

2. The model 

In our model the right-moving (left-moving) particles, represented by R (L), are never 
allowed to move towards left (right); these two groups of particles are the analogs of 
the outbound and nest-bound ants in a bi-directional traffic on the same trail. Thus, 
no U-turn is allowed. In addition to the TASEP-like hopping of the particles onto the 
neighboring vacant sites in the respective directions of motion, the R and L particles on 
nearest-neighbour sites and facing each other are allowed to exchange their positions, 
i.e., the transition RL — > LR takes place, with the probability K. This might be 
considered as a minimal model for the motion of ants on a hanging cable as shown in 
Fig. [I] When a outbound ant and a nest-bound ant face each other on the upper side 
of the cable, they slow down and, eventually, pass each other after one of them, at least 
temporarily, switches over to the lower side of the cable. Similar observations have been 
made for normal ant-trails where ants pass each other after turning by a small angle 
to avoid head-on collision fT7\ I16j. In our model, as commonly observed in most real 
ant-trails, none of the ants is allowed to overtake another moving in the same direction. 
Ants drop a chemical (generically called pheromone) on the substrate as they move 
forward ■ They can follow the trail by picking up the "smell" of the trail pheromone 
provided the pheromone evaporates sufficiently slowly with time. We now introduce 
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a third species of particles, labelled by the letter P, which are intended to capture 
the essential features of pheromone. The P particles are deposited on the lattice by 
the R and L particles when the latter hop out of a site; an existing P particle at a 
site disappears when a R or L particle arrives at the same location. The P particles 
cannot hop but can evaporate, with a probability / per unit time, independently from 
the lattice. None of the lattice sites can accommodate more than one particle at a time. 
From now onwards, we shall refer to this model as the PRL model. The state of the 
system is updated in a random-sequential manner. Because of the periodic boundary 
conditions, the densities of the R and the L particles are conserved. In contrast, the 
density of the P particles is a non-conserved variable. The distinct initial states and the 
corresponding final states for pairs of nearest-neighbor sites are shown in figO together 
with the respective transition probabilities. 

3. Results 

Suppose iV + and AL = N — N + are the total numbers of R and L particles, respectively. 
For a system of length M the corresponding densities are c± = N±/M with the total 
density c = c + + c_ = N/M. Of the N particles, a fraction = N + /N = c + /c are 
of the type R while the remaining fraction 1 — are L particles. The corresponding 
fluxes are denoted by F±. In both the limits 0=1 and = this model reduces 
to our earlier model [TH| ITH]. motivated by uni-directional ant-traffic, which is closely 
related to the bus-route models [201 EI] an d the models of pedestrian dynamics [221 E3] ■ 
One unusual feature of this PRL model is that the flux does not vanish in the dense- 
packing limit c — > 1. In fact, in the full-filling limit c = 1, the exact non- vanishing flux 
F + = Kc + C- = -Fl at c+ + c_ = c = 1 arises only from the exchange of the R and L 
particles, irrespective of the magnitudes of f, Q and q. 

In the special case Q = q = qu the hopping of the ants become independent of 
pheromone. This special case of the PRL model is identical to the AHR model |Hj 
with g_ = = k. A simple homogeneous mean-field approximation (HMFA) yields the 
estimates 



irrespective of f , for the fluxes F± at any arbitrary c. On the other hand, the exact 
expression for the flux of the AHR model, parametrized by a rapidity variable 



F± ~c± q H (l-c) + Kc: 
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Figure 3. The fundamental diagrams in the steady-state of the PRL model for several 
different values of (a) / (for <j> = 0.5) and (b) cf> (for / = 0.001). The other common 
parameters are Q = 0.75, q = 0.25, K = 0.5 and M = 1000. 



for 1/2 < qjj < 1, where a = and the unit of elementary time scale has been set 
by choosing K = 1. A comparison of the equation (with K — 1) and the exact 
result (JHJ) in Fig. H] shows that the flux in the HMFA, although an underestimate, is a 
reasonably good approximation for all qn > 1/2. Deviation from the exact values for 
qn < 1/2 indicates the presence of stronger correlations at smaller values of qn- For the 
generic case q ^ Q, the flux in the PRL model depends on the evaporation rate / of 
the P particles. In Fig. Elwe plot the fundamental diagrams for wide ranges of values of 
/ (in Fig. Of a)) and (in Fig. Ofb)), corresponding to one set of hopping probabilities. 
First, note that the data in Figs. El are consistent with the physically expected value 
of F±(c = 1) = Kc + C-, because in the dense packing limit only the exchange of the 
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Figure 4. Fundamental diagram in the special case Q = q = qn- The continuous 
curves, marked by *, + and x, are the exact results corresponding to qn = 
0.90,0.75,0.60, respectively. The corresponding HMFA results have been shown by 
the dotted, dashed and dashed-dotted curves, respectively. The points marked by *, 
+ and x have been obtained from computer simulations. 
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Figure 5. Time gap T g between successive collisions as a function of system size M 
for Q = 0.75, q = 0.25, K = 0.5, / = 0.005, c = 0.2 and cp = 0.3 

oppositely moving particles contributes to the flux. Moreover, the sharp rise of the flux 
over a narrow range of c observed in both Fig. El (a) and (b) arise from the nonmonotonic 
variation of the average speed with density, an effect which was also observed in our 
earlier model for uni-directional ant traffic [TT^I19j. This nonmonotonicity was shown to 
be a consequence of the formation of so-called loose clusters. These are regions in space 
where the particle density c\ c is larger than the average density c, but not maximal, 
i.e. c < c\ c < 1. These loose cluster dominate the behaviour at intermediate particle 
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Figure 6. Space-time plot of the PRL model for Q = 0.75, q = 0.25, / = 0.005, 
M = 4000, c = 0.2 and (a) = 0.5, K = 0.2, (b) <j> = 0.3, K = 0.2, (c) <j> = 0.3, 
K = 0.5. The red and green dots represent the right-moving and left-moving ants, 
respectively. 

densities c and small evaporation rates / where they are formed due to effectively longer- 
ranged attractive interactions introduced by the pheromones. 

In the special limits = and 0=1, over a certain regime of density (especially at 
small /), the particles are known ^Hl El to form "loose" (i.e., non-compact) clusters 
[T§] which are defined to be regions of the system with a density q c that is larger than 
the average global density, i.e., c < q c < 1. If the system evolves from a random initial 
condition at t = 0, then during coarsening of the cluster, its size R(t) at time t is given 
by R(t) ~ t 1 / 2 pm I2T] . Therefore, in the absence of encounter with oppositely moving 
particles, t±, the coarsening time for the right-moving and left-moving particles would 
grow with system size as r + ~ 2 M 2 and r_ ~ (1 — 0) 2 M 2 . 

In the PRL model with periodic boundary conditions, the oppositely moving 'loose" 
clusters "collide" against each other periodically. Let us define r g to be the time gap 
between the successive collisions of the two clusters which is the time interval between 
the end of one collision and the beginning of the next collision. It is straightforward to 
see that r g increases linearly with the system size following r g ~ M; we have verified 
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(a) (b) 

Figure 7. Magnified view of the first and second collision area respectively of Fig. 
He) 
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(a) (b) 

Figure 8. Space-time plot of the PRL model for Q = 0.50, q = 0.25, / = 0.005, 
c = 0.2, (p = 0.3, K = 1.0 and (a) M = 1000, (b) M = 4000. The red and green dots 
represent the right-moving and left-moving ants, respectively. 

this scaling relation numerically as in shown in Fig. EJ During a collision each loose 
cluster "shreds" the oppositely moving cluster; both clusters shred the other equally if 
= 1/2 (Fig. HHa)). However, for all (j) ^ 1/2, the minority cluster suffers more severe 
shredding than that suffered by the majority cluster (Fig. E(b)) because each member 
of a cluster contributes in the shredding of the oppositely moving cluster. 
Fig. Ufa) andllfb) show magnified view of the first and second collision area respectively 
in the space-time plot Fig. |Hfc) where only right moving ants are present. Fig. [7(a) and 
UJb) demonstrate that the phenomenon of shredding is such that the total number of 
clusters, their sizes as well as their individual velocities just after the collision are usually 
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Figure 9. Fundamental diagrams for the system size M = 1000 and M = 4000; 
common parameters being c = 0.2, <j> — 0.3, Q — 0.75, q = 0.25, K = 0.50 and 
/ = 0.005 




Figure 10. Average size of the cluster R plotted against time t for <fi = 1.0, and 
<f) = 0.5, both for the same total density c = 0.2; the other common parameters being 
Q = 0.75, q = 0.25, K = 0.50, / = 0.005, M = 4000. Dotted line corresponds to a 
slope of i 1 / 3 and dashed line corresponds to a slope of t 1 / 2 . 



different from those of the clusters before the collision. But, at present, we do not have 
any theory to predict the changes in these properties caused by the collision. 
In small systems the "shredded" clusters get opportunity for significant re-coarsening 
before getting shredded again in the next encounter with the oppositely moving particles. 
But, in sufficiently large systems, shredded appearance of the clusters persists. This 
is demonstrated clearly by the space-time plots for two different system sizes in 
Fig. El However, we observed practically no difference in the fundamental diagrams 
for M = 1000 and M = 4000 (see Fig. EJ). Following the methods of ref. [2H, we 
have computed R(t) starting from random initial conditions. The data corresponding 
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to = 1 are consistent with the asymptotic growth law R(t) ~ t 1 ^ 2 . In the begining 
R(t) grows as t 1//3 (dotted line in Fig. corresponds to a slope of t 1//3 ) however in the 
later stage it grows as t 1 / 2 (dashed line in Fig. El corresponds to a slope of t 1//2 ). In 
sharp contrast, for = 0.5, R(t) saturates to a much smaller value that is consistent 
with highly shredded appearance of the clusters in Fig. Efa). 

Thus, coarsening and shredding phenomena compete against each other and this 
competition determines the overall spatio-temporal pattern. Therefore, in the late stage 
of evolution, the system settles to a state where, because of alternate occurrence of 
shredding and coarsening, the typical size of the clusters varies periodically. Moreover, 
comparing Fig. Efb) and Fig.Efc), we find that, for given c and 0, increasing K leads to 
sharper speeding up of the clusters during collision so long as K is not much smaller than 
q. Both the phenomena of shredding and speeding during collisions of the oppositely 
moving loose clusters arise from the fact that, during such collisions, the dominant 
process is the exchange of positions, with probability K, of oppositely-moving ants that 
face each other. 

4. Conclusions 

The PRL model reported in this paper, is motivated by bi-directional ant traffic. In a 
spatially constrained situation, e.g., on a hanging cable (FigHJ), such a single-lane model 
is adequate, whereas otherwise a full two-lane model [24J is required. 

The main effect of the new species of particles P is that coupling of its non- 
conserving dynamics with the conserved dynamics of the L and R species of particles 
gives rise to an effective pheromone-mediated interaction between the particles of the 
same species. This pheromone-mediated interactions between the L (R) particles gives 
to a left-moving (right-moving) cluster. This tendency for "coarsening" , induced by the 
pheromones, then competes with the "shredding" process which occurs during collision 
of oppositely moving clusters in a finite system with periodic boundary conditions. 

The most surprising finding is a nontrivial and, at first sight, counter-intuitive, 
system size dependence of the spatio-temporal organization which, we believe, could be 
relevant also for other driven systems with competing aggregation and disintegration. It 
would be interesting to obtain a more quantitative description of the shredding process 
that allows to quantify this size dependence. Work in this direction is currently under 
progress. 
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